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Abstract 

The inverse problem of recovering an electrical conductivity of the 
form 7 (x) = l + ( jfe - 1)xd (xd is the characteristic function of D ) on a 
region ft C 1R 2 from boundary data is considered, where D CC ft and 
k is some positive constant. A linearization of the forward problem is 
formed and used in a least squares output method for approximately 
solving the inverse problem. Convergence results are proved and some 
numerical results presented. 
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1 Introduction 


Impedance tomography seeks to recover information about the internal elec- 
trical conductivity of an object by means of voltage and current flux measure- 
ments made on its boundary, and so provide a non-invasive, non-destructive 
imaging technique. The goal of this paper is to provide a method for the 
approximate recovery of certain types of perturbations of a constant back- 
ground conductivity. By using methods adapted to the class of conductivities 
at hand, one hopes to achieve more modest computational loads than more 
general methods (e.g., [2]) and better continuous dependence of the resulting 
estimates on the boundary measurements. 

The problem may be formulated mathematically as follows. Let J2 be a 
bounded simply connected open subset of IR^, N > 2, with C 1 2 * * * * 7 boundary 
and D an open subset of (2 with D CC 12. Define the function 7(2) as 


7 (x) 


1 x E 0 \ D 

< 

k x E D. 


Let u(x) be the solution to the elliptic boundary value problem 


L-,u = V • 7 Vu = 0 in 12 

d v u\an = 9 (1.1) 

f udS — 0, 

Ja n 

where dS denotes surface measure on 30, v is the outward unit normal vector 
to dQ y d v u = v ■ Vu, and g is a function on d£l with f an gdS = 0. Physically, 

7 represents the conductivity of the body 12, g is an applied current flux 
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density on dQ, and u denotes the electrical potential induced on 0. It is well 
known that (1.1) has a unique solution. In this setting the goal of impedance 
tomography is to recover information about the region D given the value of 
the potential u on the boundary of fl induced by the known current flux g. Of 
course one might apply several different current fluxes, measure the potential 
that each induces on dQ, and use all of this information to try to recover D. 
For notational simplicity we will consider only the case of a single boundary 
measurement, i.e., one applied current flux and measurement of the induced 
potential on dQ. It is straightforward to extend the results to multiple applied 
fluxes. The constant k will be considered known a priori. Friedman [5] has 
proven results regarding the detection and identification of the region D from 
a single boundary measurement and some results concerning the continuous 
dependence of D on the boundary data also exist (see [3]). It is known that 
D is uniquely determined if one takes all possible boundary measurements, 
that is, if one applies all possible currents fluxes and measures the induced 
voltages for each (see [7] or [8]). 

It will be necessary to restrict the domains D to lie in a certain admissible 
class. It will be assumed that this class is described by a finite number of 
parameters, so that D — D(q ) with q E Q, where Q is a compact subset of 
]R m . Certain restrictions on the map q — ► D(q) will be made later. 

We use W fc (ft) to denote the Sobolev space obtained by completing 
C°°(Q) with respect to the norm 

ll^ilw*(n) = X) l|d a< £IU J (n)- 

M<* 

The solution to (1.1) lies in W^fi) (see [6]) for sufficiently regular g. 
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In practice one does not measure the potential at every point on 
but only at finitely many locations. We will thus specify an observation 
mechanism for the potential on d£l. Let Ni, i — 1, . • • > M , be open subsets 
of dti. For each N, let fi{x) be a bounded measurable function supported in 
N{. We will then assume that the measurement of the potential on dQ, is of 
the form y = T{u) where the map T : L 2 (dCl ) -> is defined by 

H4)= \ M dS > 

JNi 

i = 1, . . . , M. Note that the restriction of a function in W 1 ^) to dU makes 
sense as an element of 2/^(c?fI) so that the integral is defined and finite if 
<f> G W 1 (fl). Since the trace operator is continuous from to L 2 (dU), 

the operator T is continuous from W 1 (Q) to lR m . 

One method for approaching the inverse problem is to use an output 
least squares method. Specifically, let di denote the potential observed by 
the electrode at Ni for some unknown region D(qo) when the current flux g 
is applied and let y® = where tt 9 solves (1.1) with D = Define a 

fit-to-data function J(q ) by 

1 M 

J(q) = ~ d i\ 2 - 

z t=i 

One may then attempt to solve the inverse problem by seeking a solution 
to the optimization problem of finding q* (E Q which minimizes J(q)- One 
of the drawbacks of this approach is that every evaluation of J(q) requires 
one to solve the forward problem (1.1), which may be costly computationally. 
Also, one has no (obvious) direct way to compute the derivatives of J(q) with 
respect to q\ availability of these derivatives is helpful in any optimization 
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problem. We will instead consider a modification of the inverse problem, 
replacing the solution u of (1.1) by u, its tangent line approximation in the k 
variable about the point k = 1. The function u may be computed explicitly, 
much more quickly than u, and using u in place of u in the definition of 
J allows us to compute directly the derivatives of J with respect to the 
parameter q, provided the map q —* D(q) is suitably restricted. 

The organization of the paper is as follows. The second section deals with 
preliminary results concerning the Neumann function for the Laplacian on a 
bounded region. The third section establishes the differentiability of the map 
k —*■ u(x ) for x € dCl for u solving (1.1) and the accuracy of the linearization 
about k = 1. The original inverse problem is replaced by an inverse problem 
involving the linearized solution u. In the fourth section a computational 
method for solving the linearized version of the inverse problem is presented 
and convergence results are proved. The final section contains numerical 
results; the algorithm for solving the linearized inverse problem is applied to 
boundary data generated by solving the linearized forward problem (for it) 
and the original equations (1.1). 


2 The Neumann Function 

The Neumann function for the Laplacian on fi is a function N(x,() defined 
on fi x f l which satisfies, for each x (E fi, the conditions 

A(W(x,0 = 

& f N(x,£)han = 


S x 

-1 

|M| 
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with the normalization /an N( x ,£) dS^ = 0. Here 8 X denotes the delta func- 
tion at x and is the Laplacian applied in the ( variable. For 0 CC IR 2 
one can verify that the Neumann function is given by 

tf(x,o = r(|s -£!) + *(*, 0 

where T(r) = — ^log(r) and a is chosen to solve 

A f cr(x,0 = 0 

d vt (T(x,£)\teBO = — — 

/ a(x,()dS i = 0. 

Jq n 

Also, if tt 0 G C 2 (Cl) satisfies 

Au 0 = / 

dv u o\dd = 9 

[ uodS = 0 

Jd n 

then uq can be represented as 

«.(*)= [ N{x,(,)g(i)iS t - [ N(x,()fU)dt (2.1) 

Jan J n 

for each x 6 fi. 

The remainder of the paper will be limited to the case in which H C IR 
and we shall next make a restriction which, for simply connected ft, represents 
no loss of generality. Specifically, we assume that ft is the unit disk, for one 
can always map the region f 2 to the unit disk conformally via a mapping 
<f>. It is straightforward to verify that the function u o <f>~ 1 defined on the 
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disk satisfies an equation of the form (1.1) with the same constant k and D 
replaced by (f>(D). The conformal map <j) also maps the unit normal vector 
field on 5ft to a normal (though not unit) vector field on the boundary of the 
disk, so that the transformed problem has Neumann data at the boundary. 
In the case in which ft is the unit disk one can verify that the Neumann 
function is given by 

W(«.0 = -^(l°g(l* - fl) + log(|* - f I)) 

where x = x/|®| 2 . 

3 A Linearization 

For a particular region D in ft and fixed Neumann data g, let u(k, x) denote 
the parameter-dependent solution to (1.1). Let Uo(x) denote the harmonic 
function on ft defined by (2.1) and define 

V>(x) = - [ Uo( t)d V( N(x,£)dS ( 

J dD 

for x 6 5ft. 

Lemma 3.1 For fixed x 6 5ft the map k — > u(k,x) is differentiable with 
respect to k at k = 1, and its derivative is given by ip(x). 

Actually, we will show something a bit stronger, namely that 
lim | u(k,x) — u 0 (x) — (k — 1)^001 < C{k — l) 2 

k — ► 1 
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for some constant C , so that the error in the tangent line approximation is 
0((k - l) 2 ). 

Before proving Lemma 3.1 we will need a couple of facts about the func- 
tion u solving (1.1). It is shown in [9], section 16 (“Diffraction Problems”) 
that u £ V^ 2 (D \ D) and u £ W 2 {D). It is also shown that the function 
u is continuous across dD and satisfies the jump condition d^u = kd x v u on 
dD, where d\,u denotes the outward unit normal derivative of u on dD as 
measured from inside D, and d^u is the outward unit normal derivative as 
measured from 0 \ D. Using these facts one can prove the following propo- 
sition. 

Proposition 3.1 For x £ Cl\D the solution u to equation (1.1) satisfies 
u(k,x) = u 0 (x) + (l-k) f u(k,t)d l/( N(x,t)dS i . (3.1) 

J dD 

For x € dD } u satisfies the integral equation 

«(M) = J aD u(k,()d,,N(x,()dS ( + (3.2) 

Proof: Since k is fixed throughout the proof, we’ll simply write u(i) for 
the function u. To prove equation (3.1), take x £ Q \ D and B t (x) a ball 
of radius e around x. Applying the divergence theorem on fi \ ( D U B e (x)) 
(valid since u is W 2 on this region) and using the fact that u(£) and JV(x,£) 
are harmonic as functions of £ on this region shows that 

0 = LM^t N M)-NM)3A0)ds ( 

J dlt 

- L, f ) - OSAO) 

- f sD (u(()dl t N(z, () - N(x,()dlu(())dS ( , 
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where the vector field v of the boundary of B € (x) points away from the ball. 
Using the properties of the Neumann function and the representation formula 
(2.1) yields 

0 = -<*>(*) - / 

- f aD m x(*,o- mx, ;m)) ds ( . 

As e goes to zero standard potential theory arguments (e.g., [4], chapter 3) 
give 

0 = -««(*) + “OO - I M()%Mx,()-N(x,()SZu({))dS ( . (3.3) 

JoD 

Note if x G dCl then B e (x) is only “half” of a disk, although in this case 
N(x,£) = 2r(x,0, 80 that (3.3) remains valid. Finally, since N is smooth 
near dD as a function of £, dl ( N — dl^N so that 

/ %0W) dS i = kf N(x,£)dlu(£) dS{ 

JdD JdD 

= k f gD d V( N( X ,()u(t)dS ( 

where the last equality follows by applying the divergence theorem on D. 
When substituted into equation (3.3), the last equality yields (3.1). 

To prove equation (3.2), rewrite the integral in equation (3.1) as 

/ u(()a,A l(x,()dS t = I Stu(()N(x,()dS ( 

JdD 5 JdD 

When this is combined with equation (3.1), we obtain, for x € 0, \ D, 

u{x) = «o(x) + ^ / dtu{ON{x,i)dS i . (3.4) 

K J dD 
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This expression remains valid for x G dD as well, since the singularity in the 
integrand is only logarithmic and u is continuous across dD. For x G dD let 
D e = D\ B f (x). An application of the divergence theorem shows that 

[ u(0d V( N(x,0dS t = f dluiONix^dSf 

JdD, JdD, 

As e tends to zero this becomes 

/ u (()d^N(x.()ds ( + hi(x) = / a>«)w(x,0rfs ( 

JdD Z JdD 

rC J dD 

or 

f dlu(()N(x i ()dS t = ^u(x) + k [ u(()d V( N(x,()dS i . 

JdD 1 JdD 

This equation, upon substitution into (3.4), yields equation (3.2). □ 

Proof of lemma 3.1: From equation (3.2) it follows that 

sup |u(&, x)| < — ~r k 3U P K*»OI)( su P \\dv( N \ Ui(ao)) + 7 TT s 1 1 pI U o I 

xedD 1 + K (effD xedD * t t S 8D 

or, with a little rearrangement, 

, ,, v. ^ 2 sup 90 |uq| 

*zbd U >X ~ 1 4- As — 2(1 — *)sup x ||a v{ lV(x f OIUj(eo)’ 

Since sup xe9D \\d V( N(x,()\\ L ^ dD ) is finite (d V( N(x,£) is actually continuous as 
a function of x and £-see [4], section 3C), sup 9I) |u| is bounded as a function 
of k for k in some neighborhood of one. In fact, rewriting equation (3.2) as 

u(k,x) - Uo(i) = J aD u(()d ul N(x,()dS t + 

and using the boundedness of u as a function of k shows that 

sup |u(A:,x) — uo(sc)| < C\k — 1|. 

x£dD 
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The constant C will depend on D and u 0 . From equation (3.1) and the 
definition of V>(x), 

u(k,x) - u 0 (x) - (k - l)i/;(x) = (1 - k) f (u-Uo)(£)d V( N(x,()dS( (3.5) 

for x 6 OU. Since sup azJ |u — u 0 | < C\k — 1| on dD, equation (3.5) shows 

sup |u(A;, x) — ii(x)| < CIA: — 1| 2 
iesn 

for x € 3f 2 and some 0, independent of x. This proves Lemma 3.1. □ 
Define the linearized solution u(x) as 

u(x) = Uo(x) + (A: — l)^»(x). 

Next we make some restrictions on the map q — ► D(q). First, we require 
that D(q) C fP for some fixed fi' CC fi and that the boundary of D(q) is a 
simple closed curve. Parameterize dD(q ) as x = x q (t),y — y q (t), 0 < t < 
with 



where L is the length of the boundary of D. We also require that for all 

qtQ, 


||x,||c*[0,2ir] < K 

l|y«llc»[0,2x] < K 

for some constant K. This follows, for example, if the curvature of dD(q ) 
is bounded. Finally, we require that if q — ► q 0 in IR ro then D(q) — ► D(q 0 ), 
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j 

where convergence of these domains means that for each e > 0 there is a j 

number 6 > 0 so that j 

ll 1 ? — X QO ||c s [0,2ir] < e > j 

l|y<i - Vqo l|c*[0,a»] < e | 

whenever — <7o| < 6 {6 is independent of go)- I 

We now reformulate the original inverse problem, using the the linearized 
solution u in place of u: 

(IDP) minimize J(q ) for q G Q 

where J(q) is defined as • 

i M 

Ew-y 1 . i 

z t=l 

but with y q = J-{u q ), where u q is the linearization of the solution to (1.1) . 

with D = D(q) and T and the are as defined in the introduction. ; 

One can prove that the problem (IDP) has a solution. The conditions on = 

the map q — > D(q) guarantee the continuity of the map 

q — ► u 0 (x ) + (k - 1) f u 0 (£)d U( N(x,£)dSt 
J dD 

as a map from IR m to L 2 (dU) (d V( N(x,£) is smooth for x 6 3f l and £ 6 dD), 
provided the Neumann data g is regular enough, for example, g € L 2 (5f2). 

The map T is continuous from L 2 (dU) to 1R M so that q — ► J(q) is continuous = 

and hence J(q) attains a minimum for some q* 6 Q. Moreover, the map 
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? -> u is differentiable if the map q -> ( x q (t),y q (t )) is differentiable, so that 
we may use an optimization scheme involving evaluation of derivatives in 
solving (IDP). 

Remark: If the Neumann data g is not sufficiently regular then the func- 
tion u Q may not be in L 3 (dU ); however, suppose that u 0 £ L 2 (S) where 
S C dCl. Then the above argument remains valid if we require that the elec- 
trode locations be subsets of S. Also, if g is smooth enough on S C dfl 
then uq will be continuous on S. In this case we may take the functions /* to 
be Sp., delta functions at points p j £ N{, so that F is simply a point-by-point 
sample. The proofs above and in section 4 remain essentially unchanged. 


4 A computational method 


In this section we will consider a computational method for solving the iden- 
tification problem (IDP). An n point quadrature rule rj n on the interval [0, L ] 
is a set of pairs of real numbers £ [0, L], for i = 1, . . . , n. We 

will say that the family of quadrature rules {r] n } = u\ n }, i = 1, . . . ,n, is 

convergent if for any <j> £ C l [0, L\, 


" rh 

/ <f>(t)dt 

i = 1 J ° 


0 


(4.1) 


as n — » oo. The convergence of the family {»/„} will be called uniform if the 
convergence in (4.1) is uniform over the set {(f> £ C ,1 [0,L]; ||<^||cqo.L] < A} 
for any constant A, i.e., the rate of convergence depends only on C l bounds 
for (f>. 
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Given a quadrature rule rj n = {i;„, u>;„} define J n (q) as 

■> M 

ru) = ± E is? - *1*. 

Z »= 1 

with y 9 = T and d{ as previously defined and u’ defined by 

u«(x) = tio(x) + (& — 1 )t— 

^ .'=i 

where £(i) is a parameterization of dD{q ) satisfying (3.6). Define the discrete 
approximation to the identification problem (IDP): 

(AIDP)" minimize J n (q) for q G Q. 


Theorem 4.1 Let {rj n } be a uniformly convergent family of quadrature rules 
and let q n be the solution to the associated problem (AIDP/ 1 . Then the se- 
quence q n contains a subsequence q nm converging to some q* 6 Q- Moreover, 
q* is a solution to (IDP). 

Proof: If D = D(q ) for some q e Q and dD is parameterized in (£i> 6) co- 
ordinates as £( t ) = (£i(i),6(t)), subject to the condition (3.6) then, written 
out in terms of the parameterization, 

V,(x) = J* uoU(t))V(N(x, m) • § dt (4.2) 

for x E dD. Since N and its derivatives are continuous as functions of x and 
£ for x 6 dD and ( € D', we can bound them uniformly for x and ( in these 
sets. Similarly uq and its derivatives are bounded on D' . Using these facts, 
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as well as D(q') C fl* for all q £ Q , one can bound the integrand and its first 
derivative (as functions of t), uniformly for q £ Q and x £ dCl. It follows 
that 

X>i»X0 (((lfa))V,W(*, «<(.)) ■ §((<„) - *(*) 

1=1 

as n -*• oo, uniformly for q £ Q, x £ so that (x) -► u q (x). The operator 
T is continuous from L 2 {dVl) to 1R M so that J n (q ) — ^ J(q) uniformly with 
respect to q. Hence if q m is a sequence in Q with q m — > q then 

= J (l)- ( 4 -3) 

Let q n be a solution to (AIDP) n . Since Q is compact some subsequence 
9 nm converges to q* £ Q. For any q E Q we have 

J n ( 9 n ) < 

Taking the limit over n and using (4.3) shows that 


so q* solves (IDP). □ 


■>(?•) < J(«) 


5 Numerical Results 

We will begin with a simple case. As mentioned, fl will be taken to be the 
unit disk in 1R 2 . The family of subdomains D(q) will consist of disks centered 
at (a, b) with radius r, so that q = ( a,b,r ). The set Q will consists of those 
points satisfying y/a 2 b 2 + r < p < 1, so that each domain is contained 
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in Q and is bounded away from dCl. It is straightforward to verify that this 
family satisfies the conditions stated in section 3. Moreover, we can directly 
compute the derivatives of u with respect to a, b, and r. 

For the applied current flux we take g = —nd T S u where d T is the counter- 
clockwise tangential derivative and S u is the dirac delta function at w € dCl. 
The harmonic function Uq with this Neumann data is known in closed form 
and for u — (1,0) is given in ( x,y ) coordinates by 

“o (*,») = jz tj— 2 - 

(1 - x) 2 + y 2 

Note that this function is not L 2 on d£l } but is smooth and bounded away 
from the point (1,0) (see remarks at the end of section 3). The electrode 
locations will be taken as pi = (cos(27ri/16), sin(27ri/16)) for i = 1,...,15 
and for the functions we simply set /j = S Pi , so T{u) is just a point-by-point 
sample of u at the points p;. 

A 10-point Gauss-Legendre integration rule was used to approximate ip(x) 
and so obtain the discrete approximation J(q) to J(q). The derivatives of 
J(q) were computed by simply differentiating the sum J(q) with respect to 
q. The function J(q) was minimized using a Levenberg-Marquardt optimiza- 
tion procedure with initial guess q = (0,0, 0.1). The simulated data for the 
linearized inverse problem was generated by evaluating u(x) at the points 
x = pi . To simulate noise in the data, independent gaussian random vari- 
ables were added at noise levels of 0%, 5%, 10% and 20%, e.g., for 10% noise 
the random variables added had a mean square value of 0.1 x ||u— Uo||ia(an) in 
the data. The results are summarized in table 1. The value of k was chosen 
as 1.1, although for the linearized problem the value of k is irrelevant. 
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Figure 1 shows the accuracy of the the linearized solution as an approx- 
imation to the true solution for As = 1.1 and k = 2.0, the region D being 
that used for figures 3 through 7. The graphs are of u — Uq and u — u 0 . 
The forward problem was solved by using Nystrom’s method (see [1]) on the 
Fredholm equation (3.2) to solve for u on dD and then representing u on 
the boundary of by means of equation (3.1). Tables 2 and 3 summarize 
the results when the optimization is applied using data from the full non- 
linear forward problem. Here the value of k is relevant; the algorithm should 
perform best for k near 1. The tables present results for k — 1.1, and k = 0.5. 

To apply this method to more general regions we need a way to describe 
D using finitely many parameters. We’ll assume dD(q ) can be described as 

x{t ) = a + r(t)cos(f) 

y(i ) = b + r(t)sm(t) 

for 0 < t < 27r, where a and b are constants and r(t) is a C 2 cubic spline 
satisfying r(i*) = i = 1, . . . ,m with tj = 2’Kijm with the endpoint condi- 
tions r(0) = r(27r) and r'(0) = r'( 27r). It can be verified that the conditions 
(3.6) are satisfied provided the qi are bounded away from zero. This algo- 
rithm was implemented as follows. First, apply the algorithm which assumes 
D is a disk to recover that disk centered at (a, b) with radius r which best 
fits the boundary data. Second, with this choice for a and b and the initial 
guess qi = r apply the optimization procedure (Levenberg-Marquardt) in the 
variables qi,. . . ,q m . A smooth penalty term, identically zero for qi > 0.01, 
was added to J(q) to ensure that the qi were bounded away from zero. This 
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term was of the form £J = i(/(x)) 2 with / taken to be 

{ c 3 (x — ci) 3 exp — c 2 /(x — ci) 3 x < Ci 

0 x > Ci 

for constants ci = 0.01, c 2 = 3.05e — 5, c 3 = 1.36e4. 

In Figures 3 through 8, m = 5, a = 0.2, b = 0.4, q 1 = q 2 = 0.032, 
q 3 = 0.056, and q 4 = q 6 = 0.04, respectively. The linearized forward problem 
was solved (with a dipole at (1,0)) and the optimization procedure applied 
to the resulting data. Figure 3 shows the region D(q) and the best fit disk. 
Figure 4 shows the region recovered by optimizing over q with (a, b ) fixed as 
the center of the best fit disk. Figure 5 shows the region D(q) and recovered 
estimate when 10% noise is added to the data. 

The optimization procedure for the linearized problem was also applied 
to data generated by solving the full non-linear forward problem. The region 
D(q) is chosen as in the previous example. Figure 6 shows the recovered 
estimate of D(q) when k = 1.1, using a single dipole at (1,0). Figure 7 
shows the same with k = 0.5 and figure 8 with k = 3.0. For k near one the 
estimate of D is good. For values of k farther from one the estimates of D 
are correct in location and approximate size, but the details of the shape are 
lost. 

Finally, one is not restricted to using a single applied current flux. An 
obvious extension of the previous results allows one to apply the linearized 
approximation for multiple current fluxes. The functional J(q) is then simply 
defined as the sum of the Ji(q) where Ji(q) is the cost functional for the 
ith current flux. For Figure 9 the region D(q) was chosen to be a disk 
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centered at (0.3, —0.3) with radius 0.05 and k = 1.5. In the first figure a 
single dipole current was applied and 50% noise added to the data from the 
forward problem. In the second figure 8 dipole currents were applied, the 
dipoles located at (cos(7rz/4), sin(7rz/4)), i = 0, ...,7, again with 50% noise 
added. 
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0.000 

* - 0.001 
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"o 

* - 0.002 
- 0.003 
- 0.004 

0 12 3 4 5 6 

Boundary of unit circle, 0 to 2ir 

Figure 2: Function u — u 0 (solid line) ?ind u — u 0 (dashed line) for k = 2.0. 
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9 

true 

Noise free 

5% noise 

10% noise 

20% noise 

a 

0.2 

0.2000 

0.2010 

0.2026 

0.2639 

b 

0.3 

0.3000 

0.2952 

0.2891 

0.2945 

r 

0.05 

0.0500 

0.0507 

0.0491 

0.0452 

iterations 


5 

5 

6 

6 



5.67 x 10- 19 

1.27 x 10“ 9 

5.46 x 10- 9 

6.94 x 10- 8 


Table 1: Summary of results for linearized inverse problem 


9 

true 

Noise free 

5% noise 

10% noise 

20% noise 

a 

0.2 

0.1999 

0.2028 

0.2258 

0.1582 

b 

0.3 

0.3000 

0.2987 

0.2899 

0.3144 

r 

0.05 

0.0512 

0.0506 

0.0503 

0.0550 

iterations 


5 

5 

6 

6 

J(q) 


7.28 x IQ" 19 

2.50 x lO" 9 

7.96 x 10- 9 

2.46 x 10" 8 


Table 2: Summary of results for full non-linear problem, k = 1.1. 
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0.12 0.18 0.24 0.30 


Figure 4: Actual region D (solid line) and recovered estimate (dashed line). 



Figure 5: Actual region D (solid line) and recovered estimate (dashed line ) 
with 10% noise added to input data. 
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